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Abstract. We performed Monte Carlo simulations in a bilayer system composed 
of two thin films, one ferromagnetic (FM) and the other antiferromagnetic (AFM). 
Two lattice structures for the films were considered: simple cubic sc and a body 
center cubic bcc. We imposed an uncompensated interfacial spin structure in 
both lattice structure; in particular we emulated a FeF2-FM system in the case of 
the bcc lattice. Our analysis focused on the incidence of the interfacial strength 
interactions between the films, Je(,, and the effect of thermal fiuctuations on the 
bias field, Heb- We first performed Monte Carlo simulations on a microscopic 
model based on classical Heisenberg spin variables. To analyze the simulation 
results we also introduced a simplified model that assumes coherent rotation of 
spins located on the same layer parallel to the interface. We found that, depending 
on the AFM film anisotropy to exchange ratio, the bias field is either controlled by 
the intrinsic pinning of a domain wall parallel to the interface or by the stability 
of the first AFM layer (quasi domain wall) near the interface. 
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1. Introduction 

Exchange bias (EB) is an ubiquitous magnetic phenomenon that usuaUy appears when 
two different magnetic media are in contact. Although EB can be observed in a large 
variety of non- homogeneous magnetic materials [1, 2], in this work we will focus on 
the case of a bilayer system composed of two films, one ferromagnetic (FM) and the 
other antiferromagnetic (AFM). 

Assuming that the Curie Temperature Tq of the ferromagnetic material is larger 
than the Neel Temperature Tpf of the antiferromagnetic one, and that the two films 
are magnetically coupled by exchange interactions, an unusual hysteresis phenomenon 
can be observed. If such a system is cooled down below T^r in the presence of an 
external applied magnetic field He f the hysteresis loops of the FM material evidences 
three important anomalies when compared with the loop of the single ferromagnetic 
film. First, a shift in the loop appears, characterized by a new center called the 
bias field Heb- This shift is due to the unidirectional anisotropy induced at the 
interface. Second, the width of the loop usually increases. Finally, the loop also loses 
its symmetry. As temperature increases, the bias field Heb goes to zero at certain 
blocking temperature T^, with Tb < TV, restoring the normal hysteresis loop of the 
isolated ferromagnet. 

Although this phenomenon was reported for the first time in 1956 [3] and despite 
the huge theoretical and experimental effort devoted to understanding its origins, there 
are still many controversial points concerning the underlying mechanisms responsible 
for such unusual hysteresis anomalies [1, 2, 4, 5, 6]. In particular, these controversies 
are in part related to the fact that EB has been observed in a great diversity of 
magnetic system, including for instance spin glasses, intrinsic inhomogeneous and 
nanoparticlc systems, as well as the bilaycred system analyzed in this paper. Beyond 
the theoretical interest, this phenomenon is also relevant because of its technological 
applications-for instance, in the design of magnetic sensor and magnetic recording 
media devices [2], among many others. 

As regards the case of a bilayered FM/AFM system, the spin structure at 
the interfacial planes is a main issue in developing the understanding of the EB 
phenomenon. In particular, AFM interfaces can be roughly classified as compensated 
or uncompensated^ depending on whether the nearest AFM plane parallel to the 
interface have zero net magnetization or not, respectively. Most of the earlier models 
that explain EB assume an uncompensated interfacial spin structure [4], even when 
this requirement is not always fulfilled in experiments. Actually, EB can be observed 
also in compensated interfaces, but in this case the existence of uncompensated 
domains has shown to be fundamental for the appearance of the hysteresis shift [7]. 
Furthermore, fully uncompensated interfaces can manifest a weaker EB field when 
compared with partially uncompensated or even compensated interfaces. In fact, 
experiments carried out by Moran et al. [8] and Nogues et al. [9] on Fe films grown over 
FeF2 single crystals cut along different orientations, showed that Heb is larger when 
the interface is compensated ((110) plane) in comparison with the uncompensated 
case ((010) plane). This effect is supposed to be associated with spin re-arrangement 
at the interface [2, 9] since a similar behavior was found when the AFM is a single 
crystals or a thin film. 

A key point for the understanding of the EB phenomenon on uncompensated 
interfaces is the effect of the variation of the exchange coupling between interface 
layers on the EB field. While it is difficult to control this quantity at the 
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experimental level, this problem can be handled easily using Monte Carlo simulations 
based on microscopic models. In addition, this methodology allows a detailed 
description of the interfacial spin structure together with the incorporation of thermal 
fluctuations, which are relevant for the stability and therefore the appearance of 
the EB phenomenon. For instance, thermal effects are necessary to explain the 
widening of the hysteresis loop close to the blocking temperature [10, 11, 12]. In 
this sense, numerical studies at the micromagnetic [13, 14, 15] and Monte Carlo 
simulations levels [16, 17, 18, 19, 20, 21, 22, 23] have proved to be very useful tools 
for modeling realistic systems. On the other hand, the continuous approximation 
assumed in micromagnetic based model breaks down in highly anisotropic materials 
like FeF2 antiferromagnetic compounds. Discretization could give rise to different 
energy barriers with the consequent thermal activated effects [24]. Hence, atomic 
scale based models could be crucial for getting an appropriated description of the 
magnetic properties. 

In this paper we analyzed the EB phenomenon in a FM-AFM bilayer system 
with an uncompensated interface. In section 2 we summarize the existing theoretical 
background, discussing the phenomenology of EB system in the frame of two of the 
most relevant models. In section 3 we introduce a microscopic model for the bilayered 
system, describe the simulation protocol and show our numerical results. In order to 
interpret the results of the previous section we introduce in section 4 a generalization 
of Meiklejohn-Bean model, which allowed us to analyze the role of the AFM layers in 
the EB phenomenon. In section 5 we summarize and discuss the results. 

2. Theoretical Background 

In order to analyze the role of the strength of the interface exchange interaction J^b in 
the behavior of bias field iJeb, let us discuss first the following question: what happens 
with the order of the AFM as we invert the orientation of the FM magnetization 
by applying an opposite magnetic field hi We assume that the system has already 
reached thermal equilibrium at certain temperature T below Tb, in such a way that, 
if Jeb were zero, both films would have achieved an ordered state. Since Tc > T^, 
we assume \Jf\ > \Ja\ where Jp and J a are the exchange interactions of the FM and 
AFM, respectively. If J^b is small enough {J^b << Ja) the spins in the AFM will 
remain almost insensitive to the rotation of the global magnetization of the FM film. 
In this case the Meiklejohn-Bean model [25] predicts a linear dependency of the bias 
field Heb on the value of Jeb'- 



where L pm is the thickness of the FM film. 

At the other extreme, when Jeb is large enough, the rotation of the magnetization 
would induce the creation of a domain wall (parallel to the interface) in the AFM films, 
at least for small enough values of Ka- Once a perfect domain wall has been formed, 
any increase of Jeb will not alter the value Heb- This phenomenology is captured by 
the model of Mauri et al. [26] (from now on the MSBK-model) when the anisotropy 
of the FM film is negligible. This model predicts an initial increase of Heb with Jeb 
for small values of Jeb, followed by a saturation for large enough values of Jeb at 
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where is the anisotropy constant of the AFM and w is constant depending on the 
lattice structure. The previous results suggest a monotonous behavior of the bias field 
when the anisotropy of the FM film is negligible, with a linear dependence of i?eb with 
Jeb for small values of Jeb and a saturation for large values of it. As we will show 
in the next section, such scenario can change substantially depending on the effective 
anisotropy of the AFM film. 

3. Microscopic model and Numerical Simulations 

3.1. The microscopic model 

We considered a FM film mounted over an AFM film as depicted in figure la. The films 
are magnetically coupled to each other by exchange interactions and the structure of 
both films is either hcc or sc, assuming a perfect match across the FM/AFM interface. 
The system is ruled by the following Hamiltonian, 

H = -Jp Sp- Sp, ~ Kf ('^--)' 

<f,f' >eFM feFM 

- Y JAF{r^r')Sr-Sp,~KA ^ {Slf 
<r,r'>eAFM feAFM 

-Jeb S?-Sp-hYsl. (3) 

<r,r'>eFM/AFM r 

where Sp is a classical Heisenberg spin {\Sp\ = 1) located at the node r of the lattice. 
< f,f' > denotes a sum over nearest-neighbors pairs of spins, Jp > is the exchange 
constant of the FM and JAF{r — r') is the strength of the AFM exchange interactions 
which explicitly depends on the vector r — r' . This dependency of Jaf on f—f' is 
introduced in order to set an uncompensated interface at the AFM. For the bee lattice 
we set Jaf = —Ja with Ja > for all pairs (r, r'), while for the sc lattice we set 
Jaf = J A if ('^j belong to the same plane parallel to the interface and Jaf = — Ja 
otherwise (see figure lb). Jeb > represents the exchange coupling between the FM 
and the AFM interface layers of the films, Kp is the anisotropy constant of the FM, 
Ka is the AFM anisotropy and h is an external homogeneous magnetic field oriented 
along the y direction. We assumed that: 

(i) Kp < in order to ensure the FM anisotropy term tends to align the spins on 
the plane of the film, mimicking the dipolar shape anisotropy, as usual [23, 24]; 

(ii) Ka > in order to introduce an uniaxial anisotropy along the y direction in the 
AFM material [18], 

We carried out Monte Carlo simulations using Metropolis algorithm and 
Hamiltonian (3). In our simulations and Ly are the lateral dimensions of the 
films, in units of the lattice parameter, and L^a and L^f are the thicknesses of the 
FM and AFM films, respectively, measured in the same units. The total number 
of spins is then N = rj Ly {L^a + L^f) where 77 = 1, (2) for the sc, (bcc) lattice. 
Periodic boundary conditions were imposed in the plane of the film while we used 
open boundary conditions in the perpendicular direction to the film. For each point 
in the magnetization curve presented in this work, we took 10^ Monte Carlo Steps per 
site (MCS) to thermalize the system and then the same number of MCS to calculate 
the temporal averages of the magnetization. We follow the ideas used in Refs. [27, 28], 
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Figure 1. a) Scheme of the bilayer system including the reference frame used 
throughout this paper, b) Schematic picture of the system modeled by the 
Hamiltonian (3) in the sc lattice case. Here we show ground state configuration 
with the corresponding interactions. 

where at each spin actuahzation the direction of the spin is adjusted in a cone in such 
a way to maintain an acceptation rate close to 0.46. This is an approximation to a 
Landau-Lifshit-Gilbert Langevin dynamics in the high damping limit [29]. We set the 
following dimensions for the system, = Ly ~ 40 and L^a — Lzf = 12. and fix 
the following parameters: Jp = 9.56J, Ja = —J and Kp = — 0.5J, where J is an 
arbitrary parameter that sets the energy units. Jeb varies in the interval [0, Jp] while 
Ka can take arbitrary values. With these parameters we can emulate a FeF2-FM 
system in the bcc lattice by choosing Ka = 1-77 J [18]. Since we are interested in 
the high AFM anisotropy to exchange ratio regimen, which implies small domain wall 
width, the thickness of the AFM we chose is enough to support an AFM domain wall. 
On the other hand, it is known in this model [18, 30] that for such sizes both the AFM 
and the FM films reach an ordered state. 

3.2. Results 

In Fig. 2 we present the bias field HpB (open circles) and coercivity He (squares) 
obtained from Monte Carlo simulations as function of the interfacial interaction 
strength J^b for the two considered lattice structures and for fixed values of 
temperature and AFM anisotropy. When the interfacial exchange coupling j^b = 
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Figure 2. Bias field Heb (circles) and cocrcivity He (squares) vs. at 
T I J = 0.5 and and Kj^/Jj^ = 1.77. (a) sc lattice, (b) bcc lattice. 




Figure 3. Hysteresis loops of several atomic layers of the bilayer corresponding to 
the FM planes (top) and AFM planes (bottom) for jet, = j^"^ and T/Ja = 0.5. 
Left panels ((a) and (c)) correspond to the sc lattice and right panels ((b) and 
(d)) to the bcc lattice. See text for details. 



Jeb/ Ja is weak, Heb shows, for both lattice structures, a hnear dependence, indicating 
that the AFM spins located near the interface are fixed, and the FM film reverses its 
magnetization by coherent rotation[16]. As jeb increases, the bias field reaches a 
maximum value at j^"^ and then abruptly drops to an almost constant value. Notice 
that the drop is larger for the sc lattice than for the bcc one. As it will be shown later, 
such effect is due to a reduction of the effective anisotropy of the AFM layer in the 
bcc case. 

In Fig. 3 we show the hysteresis loops of several planes of the FM and AFM films. 
These loops were obtained at jeb ~ jeb"^^' .i^^t before the drastic drops observed for 
Heb in Fig. 2, where the exchange bias effect is more pronounced and the cycles are 
still almost reversible. 

In Figs. 3a and 3b we present the magnetization in the interfacial (IL) and outer 
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Figure 4. Hysteresis loops of several atomic layers of the bilayer corresponding 
to the FM planes (top) and AFM planes (bottom) for j^i, = 2 and T/Ja = 0.5. 
Left panels ((a) and (c)) correspond to the sc lattice and right panels ((b) and 
(d)) to the bcc lattice. 



(OL) atomic layers of the FM film. These results show that, in the two lattices, the 
spins rotate almost coherently. In Figs. 3c and 3d we show the loops of the four AFM 
layers nearest to the interface (Ln stands for the n-th atomic layer). Comparing the 
behavior in both structures we see that the sc lattice is more flexible than the bcc 
inside the FM, but more rigid inside the AFM, because the effective anisotropy in 
the sc is larger. In particular, in the AFM film of the sc (Fig. 3c) only the first layer 
feels the effect of the FM film. In the bcc (Fig. 3d) we clearly see the formation of a 
quasi-domain wall. 

In Fig. 4 wc plot the same quantities as in Fig. 3 for a value of jeb above the 
peak, where the bias field has already diminished abruptly. Unlike the previous case 
(Fig. 3), here the AFM layers show hysteresis behavior for both the sc and the bcc 
lattices (Figs. 4c and 4d respectively). This indicates that the drop in the bias field 
is associated with the onset of irreversible changes in the magnetic dynamics. As 
occurred below the peak (Fig. 3) the changes in the AFM are constrained to the first 
planes near the interface. It is worth stressing that now the hysteresis phenomenon 
also appears in the AFM layers, as evidenced in the behavior of the coercivity in Fig. 2. 

Next we analyzed the temperature dependence of the overall magnetic behavior. 
In Fig. 5 we present the bias field Heb and the coercivity He as a function of 
temperature for two values of the interfacial interaction strength: jefc = 0.5 (Figs. 5a 
and 5b) and 1.0 (Figs. 5c and 5d). The left panels correspond to the sc lattice and the 
right ones to the bcc lattice. The dotted lines represent the staggered magnetization 
of the AFM at zero external magnetic field, normalized with respect to the value of 
the bias field at the lowest temperature. In Figs. 5b, 5c and 5d we observe that the 
system presents a blocking temperature Tb separating two phases each with different 
magnetic behavior. At low temperature the system is characterized by the presence 
of exchange bias and almost zero coercivity. On the other hand, for Tb < T < Tn 
the bias disappears and the He increases and further decays following the behavior 
of the normalized staggered magnetization. A complete different behavior is observed 
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Figure 5. Heb ^'Hd He vs T j J for Kj\/ Jj\ = 1.77 and two interfacial exchange 
interactions: top jet, = 0.5 ((a) and (b)) and bottom j,,;, = 1 ((c) and (d)). Left sc 
((a) and (c)) and right bcc ((b) and (d)) lattices. The dotted lines represent the 
staggered magnetization of the AFM at zero external magnetic field, normalized 
to the value of the bias field at the lowest temperature. 



0.08 I , ^ , ^ , ^ , ^ , ^ r 




Figure 6. Heb vs j^i, for T/Tc = 0.14 and different values of K' = 
Open symbols: bcc lattice. Filled symbols: sc lattice, (a) K' = 
K' = 0.4425. 
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in Fig. 5a, where the ordered phase coincides with the bias regime and therefore 
Tg ~ Tpf. In this case we do not observe any trace of coercivity in the simulations. 
Note that for both structures, sc and bcc, the blocking temperature decreases with 
the interfacial interaction strength, indicating that the energy barrier for de-pinning 
the partial domain wall decreases as the wall approaches to a 180° domain wall. 

Finally we explored the effect of the lattice structure on the bias field. The main 
difference between both lattice structures is the number of nearest neighbors belonging 
to adjacent layers of any site in the AFM film, which is four times larger in the bcc than 
in the sc structure. Hence, one would expect the effective anisotropy to be reduced by 
a factor of 4 in the bcc lattice respect to the sc one, while the bias field is expected to 
be four times larger in the bcc than in the sc. To check this hypothesis we calculated 
the bias field as a function of jet in both lattices for the same value of K' = Ka/uJa: 
with a = 4 for the bcc lattice and a = 1 for the sc lattice. In order to compare the 
results, one has to take into account that the Curie temperature is different for both 
lattice structures. Hence, both calculations were carried out keeping T/Tc constant. 
In Fig. 6 we plot Heb/a as a function of jeb for high and low values of K' . For large 
enough values of the anisotropy the previous conjecture is verified, namely, the only 
effect of changing the lattice structure is a rescaling of the bias field and the effective 
anisotropy. For small values of the anisotropy such scaling is observed as long as no 
hysteresis effects appear, namely, for small enough values of the coupling j^b- For large 
values of jeb the bias field exhibits only a small drop and it saturates at a constant 
value in the bcc lattice, but it drops to zero in the sc case. We observed that such large 
drop is due to the depinning of the quasi-domain wall, i.e. to a complete reversion 
of the staggered magnetization at the AFM film. This effect is not observed in the 
bcc lattice (at least for the range of parameter values analyzed here). It is due to 
a reduction in the in plane magnetization component at the AFM layers, associated 
with a canting of the spins which enhance the pinning of the wall. 

4. Layered model 

As we have seen in the previous section, the behavior of the bias field is strongly 
determined by the magnetization dynamics of the atomic layers close to the interface. 
Moreover, we observed that, for reasonably large values of the anisotropy the spins in 
each layer rotate almost coherently under the application of an external field parallel 
to the interface. On the basis of these observations, we introduced a generalization 
of Meiklejohn-Bean [25] model that explicitly includes the contribution of the AFM 
layers close to the interface. We consider that only the n layers of the AFM film 
closest to the interface are free to move, while the rest of the AFM layers keep the 
equilibrium antiferromagnetic configuration of the bulk at temperature T. Let 5* and 
Sj be the average magnetization per layer per unit area at the FM and the AFM j-th 
layers respectively. S and Uj (j = 1, . . . , n) are assumed to be unit vectors parallel to 
the interface. The magnetization per unit area of the whole FM film is then given by 
LpAiS (with LpM the FM film thickness), since we are assuming a coherent rotation 
of the whole FM film. The n-th layer is the closest one to the interface. We assume 
that the applied field H is parallel to the interface and only interacts with the FM 
film. This approximation is valid as far as the applied field is small enough compared 
with the molecular field of the AFM. Finally, we consider the anisotropy of the AFM 
to be much larger than the FM one, so the latter can be neglected. The Hamiltonian 
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of the model is then given by 

n n— 1 

Un = - Ka (""if + (-l)"a^Aao(r)t7f + aJA ^ 

i=l i=l 
-JebSn.S-H'.S, (4) 

where H' = LfmH, a = 4 (a = 1) for the bcc (sc) lattice and (Jo{T) is the average 
equilibrium magnetization per unit area of one layer in the AFM bulk, assumed to 
be parallel to the easy axis. The (—1)" factor in the second term of Eq.(4) ensures 
the correct equilibrium configuration at zero temperature and magnetic field with 
the n — 1 AFM spin aligned with the FM spin. The model is then equivalent to a 
n + 1-spin chain, where the first spin in the chain is subjected to a local effective field 
produced by the ordering in the AFM bulk, while the spin located at the end of the 
chain (5) represents the FM film which interacts with an external magnetic field and 
is ferromagnetically coupled to the n*^ AFM spin. 

At T = the sublattice magnetization in the bulk is saturated, so we have 
(To(r) = 1. In a first approximation we can consider the simplest case of only one 
AFM layer free to move n = 1 (see fig. 7a), which is enough to illustrate the general 
mechanism. The energy is then given by 

E = -Ka {cr^f - aJA cr^ - aJ^b a.S-H'.S, (5) 

where a = ffi. The FM and AFM spins can be expressed in term of the angles and 9 
respect to the easy axis direction y of the AFM, in our case the field cooling direction 
(Fig. 7b). Then 

E — — Ka cos^ (j) — aJA cos — aJ^b coa{9 — cj)) (6) 
-H' cos(6l-7), 

where the angle 7 gives the applied field direction (Fig. 7b). From now on, we will 
consider the applied field parallel to the easy axis direction (7 = tt). In order to 
obtain the hysteresis loops and the bias field, the model is analyzed using standard 
procedures (see e.g. Ref.[31]). First, we equal to zero the partial derivatives dgE and 
dfj,E in order to obtain the critical points: 

= aJeb sin(6l -(/))- H' sm{e) (7) 
= — aJeb sin{9 — 4>) + uJa sin (j) + Ka sin(20) 

and then we look at the stability criteria, dgge d^^e — de^e^ > and dgge > to decide 
whether is a minimum or not. It turns out that, 

< aJeb cos{9 -(f))- H' cos 9 (8) 
< aJeb cos{9 — (f)) [aJA cos + 2Ka cos 2(j)] 

— H' cos9[aJeb cos{9 — 0) + uJa cos + 2Ka cos 20] 

For = we recover to the Meiklejohn-Bean model (see Ref. [[6]] and references 
therein) and the bias field is given by. 



HU = aU\ll-[^] , (9) 



11 



(a) 



a 



AFM 

bulk/ ■ 

A 

^ Ja Jen 



s 



H 




Figure 7. (a) Scheme of the model for n = 1 (Eq.(5)). (b) Angles representing 
the state of the system. 



provided that Jet < Ka- In this range of J^b the coercivity field is zero. For 
Jeb/KA > 1 the bias field drops to zero, while the coercivity jumps to a finite value 
(See Fig. 8), due to the complete reversal of all the spins in the AFM film. In the limit 
Ka = oo Eq.(9) predicts a linear increase of 7?^^ with J^t (See Fig. 8). This case sets 
an upper limit for the bias field of any model with uncompensated interface. 

For Ka = {Ja 0) the coercivity is always zero and the bias field is given by 
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This equation is valid for any value of Jef,, showing a saturation at K'^^ — uJa for 
large values of it (See Fig. 8). Eq.(lO) becomes equivalent to the MSBK-model bias 
field with zero anisotropy at the FM film, if we replace J a by the partial domain wall 
energy, namely J a — > 2\/KaJa- In the general case when Ka the coercivity is 
non zero and the problem has to be treated numerically. 

To understand the general behavior of the bias field as a function of Jeb let us 
first analyze the structure of the energy landscape given by Eq.(6) in the absence of 
external magnetic fields. Suppose that the system was cooled under the presence of 
an external field Hcf pointing to the positive y direction. Then, the energy has an 
absolute minimum, corresponding to both magnetic variables S and a pointing to the 
positive y direction. We denote this minimum by (t,t)- If the anisotropy is weak, 
Ka < Ja/'2, this minimum is unique. When Ka > Ja/'^ a second (local) minimum 
appears corresponding to both variables 5* and a pointing to the negative y direction. 
We denote this minimum by (i,i). If Ka ^ J a the energy difference between both 
minima is AE k,2Ja- 

Let us consider now the descending branch of an hysteresis cycle, that is, we 
saturate the sample with an external field pointing to the positive y direction and 
decrease the field at regular steps until the sample is saturated in the opposite 
direction. The effect of the inverse applied field on the magnetic configuration depends 
on the relative strength of j^b = Jeb/ J a- If Jeb ^ the FM layer aligns with the field 
when h = H' / Ja ~ jeb but the AFM layer still points up, that is, the lower minimum 
(t,t) changes its configuration to (t,i)- Therefore, heb ~ jeb- When jeb ~ 1 (and 
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Figure 8. Bias field H'^^ as function of interfacial exchange constant J^i, for 
different values of J a and Ka- (1) Ja = and ICa > 1 (Eq.(9)); (11) Ja = and 
Ka = 1 (Eq.(9)); (111) Ja = 1 and = (Eq.(lO)); (IV) Ja = 1 and Ka = 1. 



therefore heb ^ 1), the second minimum corresponding to the configuration 
becomes absolute. As jeb further increases, the configuration (f, i) remains as a local 
minimum, until above certain value of jet it losses stability. Hereafter we will consider 
a = 1 {sc lattice) for simplicity. 

The typical behavior of the bias field for finite values of Kj^ and is illustrated 
in Fig. 8. For low values of Jeb the bias field shows a monotonous behavior, taking 
values between those given by Eqs.(9) and (10). At this regime, the local minimum 
(t, i) of the energy remains stable and the AFM layer forms a reversible quasi-domain 
wall, without inversion of its magnetization. Above some maximum value J^"^, the 
local minimum losses stability giving rise to an irreversible inversion of the AFM 
layer magnetization and the system exhibits finite coercivity and a sudden drop in 
the bias field. However, at odds with the Ja = case, the bias field drops to a 
finite value, after which it increases again monotonously with J^b (in agreement with 
the simulation results of the previous section), due to the competition between the 
anisotropy and interaction of the AFM layer with the AFM bulk magnetization. For 
large enough values of J^b the bias field saturates into a smaller value than the Ka = 
case (i?gj ^ Ja)- As Ka increases both the drop in the bias field, as well as the value 
of Jeb where it happens increase. 

Next we compared the predictions of the model with the Monte Carlo results. 
In Fig. 9 we illustrate the typical behavior for large values of the anisotropy. We 
see that, as temperature fluctuations decrease, the maximum in the bias field as 
well as the value of jeb where it occurs increases, due to thermal activation. Of 
course, this depends on the time scales involved in the hysteresis loop, i.e., on the 
rate of variation of the field. If the rate of variation of the field is kept constant, the 
Monte Carlo results systematically approach the behavior predicted by the model as 
the temperature decreases, because the characteristic activation time systematically 
increases. 

The range of anisotropy values for which the present approximation applies 
can be estimated as Ka/ Ja > f since it is known in this range the domain wall 
width is equal to one lattice parameter [32]. When the anisotropy decreases, the 



13 



— Model T=0 



X 



0.1 



CQ 



0.15- 



0.05- 








m 



Figure 9. Bias field Heb jeb ior different values of the temperature in the 
sc lattice the anisotropy value Ka_/Ja = 1.77. 



domain wall thickness increases and more layers have to be considered for a proper 
description. For small enough values of the anisotropy a smooth domain wall is 
expected, so the behavior of the system should be well described by the MSBK model. 
The crossover to the regime of the MSBK model behavior can be estimated as the 
point where the energy of the domain wall equals the exchange energy of the AFM. 
namely 2\/KaJa — Ja, which corresponds to K^/Ja = 0.25. This is illustrated 
in Fig. 10, where we compare the maximum bias field h]^g^ ~ H^'^^ /Ja and the 
minimum after the drop /i^]^ = / Ja (see Fig. 8) with the bias field predicted by 
MSBK model hsB — '2,\/Ka/ Ja- The vertical dotted lines divide the graph in three 
regions of qualitatively different behavior. The region of validity of the present model 
(Ka/Ja > §) is marked as III. In this region a quasi-domain wall forms and, unlike 
for the continuous approximation where the internal domain wall spins change their 
orientation in a reversible way, now these spins can have an irreversible or hysteretical 
behavior, like when defects are present in the AFM [24]. 

In region I the continuous approach assumed in the MSBK-model is valid. In 
region II neither the present model nor MSBK model are expected to be valid, since 
the micromagnetic approach fails because the magnetization profile is no smooth in 
the atomic scale, but more than one interfacial plane is involved in the magnetization 
process at the interface. Moreover, we have seen from the Monte Carlo simulations 
that in this region lattice structure effects can be very important. At the crossover 
point Ka/Ja = 0.25 we sec that h^g^ = , i.e., hysteresis disappears as expected. 
It is worth noting that size effects in the AFM become relevant only in regions I and 
II, in particular when the thickness of the AFM is comparable to the domain wall size. 

Let us analyze thermal effects in the bias field when the AFM domain wall is 
pinned i.e., it is not able to propagate in the bulk of the AFM material. Suppose that 
the rate of variation of the field is small enough so that the system can be assumed 
at thermodynamical equilibrium at every step of the loop. Then, the equilibrium 
behavior can be obtained by computing the partition function 
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Figure 10. Reduced bias fields /i^^^ and h'^^g as a function of the reduced 
anisotropy Ka/Ja- The dashed line is given by the MSBK-modcl (hEB = 
2^Ka/Ja) 
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dOe^'^-^, (11) 







where we have taken uJa = 1, = ^/ksT, (pi and 9 are the angles with respect to the 
y axis of the i — th AFM and the FM spins respectively and uj = H' + J an [J = a Jeb)- 
We assumed that the bulk AFM magnetization per layer is given by the mean field 
approximation [10, 33], namely 

ao(r) ^ /:(z/3ao(T)), (12) 

where C{x) is the Langevin function and z is the number of nearest neighbors which 
depends of the lattice structure z = 6, (8) for the sc, (bcc). 

The last integral in Eq. (11) can be easily solved obtaining the general expression 
(aside from an irrelevant multiplicative factor): 

Z„ = / d(j)i--- / T.?=i':os^<l>i-(-'^V'^o(T)cosc/>i-J27=i cosi!t>i-<pi+i)) 

Jo Jo 

X /o(/3w(0„)), (13) 

where /i/(x) is the modified Bessel function and w(^) = H''^ + J'^ + 2 J H'cos (p. 
The average magnetization in FM layer can be obtained as = {cos 9) = 
and the magnetization at the j — th AFM layer, mf^ = {cos can be computed in 
a similar way. Solving numerically the previous equations as function of the applied 
field and temperature we obtained the dependency of the bias field on temperature. 
We considered the cases n = 1 and n = 2. No qualitative differences were observed. 
We present here the results for n = 1, which are adequate for illustrating the general 
behavior. 

In Fig. 11 we compare the equilibrium reduced bias field Heb sls function of jeb at 
low temperatures (full lines) with the zero temperature curves obtained from Eq.(5) 
(dotted lines) for several anisotropy values. One can see that for low intcrfacial 
interaction strength jeb <C 1 the temperature has little effect on the bias field. In 
both cases an increase in the anisotropy enlarge the range of the linear behavior 
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Figure 11. Reduced bias field Heb fs. jeb for difli'crent values of Ka/Ja- Full 
lines: equilibrium eurvcs for n = 1 at T/Ja = 0.1. Dotted lines: T/Ja = 0. 

expected in the strong anisotropy limit (see fig. 8). This can be easily understood if 
we recall that in this regime the system behaves reversibly even at zero temperature. 
In other words, in both cases the behavior of the system is governed by the absolute 
minimum of the energy, so the relation jeh ~ hgh still holds, no matter the value of 
the anisotropy is. 

The main difference appears for high values of jet- First of all, the drop in Heb 
observed in the T = curves is absent in the thermalized curves, since of course at 
equilibrium there is no coercivity. Second, the bias field Heb saturates to the value 
fiEB ~ 1 as jeb increases (jet > 1) independently of the anisotropy, contrasting with 
the T = curves where the maximum value of Heb increases with the anisotropy. 
When jeb S> 1 the applied field changes the relative depth of the two energy minima. 
When ft, ~ 1 the two minima have the same energy and the magnetization at the FM 
layer inverts = 0, independently of the anisotropy. Therefore, heb ~ 1, i-e. the 
bias field reaches the saturation value observed in fig. 11. On the other hand, the bias 
field at zero temperature continuously grows with the anisotropy due to the fact that 
the energy barriers between the minima increase with the anisotropy. 

It is worth remarking that, even at equilibrium, the bias field exhibits a maximum 
at jeb 1 for large values of the anisotropy. 

5. Discussion 

We found that, in fully uncompensated interfaces, the bias field displays a non 
monotonous dependence on the interfacial interaction strength. Depending on the 
temperature and on the anisotropy to exchange ratio Ka/Ja of the AFM, Heb 
can present a peak as a function of Jeb- In particular, the peak is observed at low 
temperatures and high enough ratios Ka/Ja- When it is present, the peak position 
moves toward lower values of Jeb as the temperature is increased, while below a certain 
temperature (low compared with the blocking temperature) the peak disappears. The 
peak is associated with the onset of coercivity, i.e., with the appearance of hysteresis 
for large values of Jeb- 
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When Ka/Ja > | (region III in Fig. 10), the behavior of the bias field is 
completely determined by the dynamics of the interfacial AFM layer. For low values 
of Jeb the interfacial layer rotates coherently forming a quasi-domain wall that changes 
reversibly with the applied field. In this regime the bias field increases almost linearly 
with Jeb and thermal effects are negligible. Above a certain critical value of Jeb the 
quasi domain wall loses stability and the magnetization of the interfacial AFM layer 
changes irreversibly. In other words, the bias field is controlled by the stability of the 
interfacial layer. This scenario, supported by both the Monte Carlo simulations and 
the simple layered model introduced here, explains why the bias field can be drastically 
reduced by increasing the interfacial interaction strength (Fig. 2). Also in this regime 
{Ka/Ja > |), the behavior of the bias field is independent of the lattice structure. 
In other words, a change in the crystalline structure is just equivalent to a rescaling 
of the effective anisotropy of the AFM. 

When Ka/Ja < |, the system can still exhibit hysteresis and a peak in the 
bias field (region II in Fig. 10), but the the width of the domain wall increases as 
Ka/Ja decreases. In this case bias, field is controlled by the intrinsic pinning due to 
the anisotropy. Namely, for large values of Jet the bias field reduces because of the 
depinning of this domain wall, which depends strongly on the lattice structure. In 
particular, preliminary results showed that the pinning is stronger in the bcc than the 
sc lattice, due to canting effects in the AFM layers. A detailed study of such effect is 
underway and will be published elsewhere. 

In both regimes (II and III) the maximum bias field is smaller than the value 
predicted by MSBK-model. These results offer certain insights about experimental 
findings in FeF2 systems[8, 9] {Ka/Ja > |), where in a fully uncompensated interface 
the bias field is much lower than expected. In particular, it becomes noticeable at 
very low temperatures. According to our results, if the interfacial strength interaction 
is strong the bias field becomes no null only at very low temperatures compared with 
the Neel temperature of the antiferroagnet. 

Summarizing, depending on the anisotropy to exchange ratio Ka/Ja the bias 
field is controlled either by the intrinsic pinning of an extended domain wall parallel 
to the interface (low anisotropy regime) or by the stability of the first AFM interfacial 
plane near the interface (sharp domain wall limit). 
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